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ABSTRACT 

We present a series of numerical simulations that explore how the ‘X-factor’, Xqo - the 
conversion factor between the observed integrated CO emission and the column density 
of molecular hydrogen - varies with the environmental conditions in which a molecular 
cloud is placed. Our investigation is centred around two environmental conditions in 
particular: the cosmic ray ionisation rate (CRIR) and the strength of the interstellar 
radiation field (ISRF). Since both these properties of the interstellar medium have 
their origins in massive stars, we make the assumption in this paper that both the 
strength of the ISRF and the CRIR scale linearly with the local star formation rate 
(SFR). The cloud modelling in this study first involves running numerical simulations 
that capture the cloud dynamics, as well as the time-dependent chemistry, and ISM 
heating and cooling. These simulations are then post-processed with a line radiative 
transfer code to create synthetic ^^CO (1-0) emission maps from which Xqo can be 
calculated. We find that for 10^ Mq virialised clouds with mean density 100 cm“^, 
Xqo is only weakly dependent on the local SFR, varying by a factor of a few over 
two orders of magnitude in SFR. In contrast, we find that for similar clouds but with 
masses of 10^ Mq, the X-factor will vary by an order of magnitude over the same 
range in SFR, implying that extra-galactic star formation laws should be viewed with 
caution. However, for denser (10^cm“^), super-virial clouds such as those found at 
the centre of the Milky Way, the X-factor is once again independent of the local SFR. 
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1 INTRODUCTION 


In both the Milky Way, and in other local galaxies, star for¬ 
mation is known to take place within large clouds of molec¬ 
ular gas, the so-called giant molecular clouds (GMCs). Un¬ 
derstanding the properties of these GMCs is important for 
the light it sheds on the process of star formation. Unfortu¬ 
nately, the main chemical constituent of the gas in GMCs, 
molecular hydrogen (H 2 ), is extremely difficult to observe in 
situ, owing to the fact that the characteristic temperature of 
the gas in typical GMCs (10-20 K; Bergin Tafalla[2007 ) is 


much smaller than the temperature required to excite even 
the lowest rotational transition of the H 2 molecule. There¬ 
fore, most observational studies of GMCs rely on carbon 
monoxide (CO) as a proxy for molecular hydrogen. Obser¬ 
vations of nearby GMCs find a surprisingly tight correlation 
between the integrated intensity coming from the clouds in 
the J = 1 ^ 0 rotational transition line of ^^CO and the 
molecular hydrogen column density of the clouds (see e.g. 


Dickman||1978| |Sanders et al. 1984 

ISolomon et al.||1987| 

Strong & Mattox||19961 |Dame et al. 

200 1 | or the recent re- 


view by Bolatto, Woffire & Leroy 

2013 

), and yield a conver- 

sion factor between CO intensit} 

^ and H 2 column density 

given approximately by (Dame et al.|2001) 


Xco,gai = = 2 X K-^ km”^ s, (1) 


where VFco is the velocity-integrated intensity of the CO 
J = 1 ^ 0 emission line, averaged over the projected area 
of the CMC, and iVH 2 is the mean H 2 column density of the 
CMC, averaged over the same area. This latter quantity can¬ 
not be directly determined from observations of H 2 emission, 
but can be inferred from measurements of the total column 
density if the Hi column density has also been mapped. The 
total column density of gas in the cloud can itself be deter¬ 
mined for nearby clouds using extinction mapping (see e.g. 
Kainulainen et aL]|2009| [Froebrich fc Rowles||2010l |Pineda| 

et al.|| 2010 | ) or measurements of the diffuse 7 -ray flux pro- 


^ Unless otherwise noted, when we refer simply to CO from this 
point on, we mean ^^CO. 


© 0000 RAS 
























2 


duced by interactions between high energy cosmic rays and 


atomic hydrogen, atomic helium and H2 (e.g. Digel et ah 


The fact that the values for Xco derived for nearby 
clouds show very little variation has led to its widespread use 
as a general CO-to-H2 conversion factor (commonly referred 
to as the ‘X-factor’), even in environments very different 
from the Galactic interstellar medium (ISM). It is therefore 
important to understand when the adoption of a constant 
value for Xco is appropriate, and when it may be seriously 
misleading. 

In clouds with a fixed H2 column density, any variation 
in Xco that occurs from cloud to cloud must be due solely 
to a variation in Wco- Numerical simulations show that this 
is sensitive to three main parameters: the velocity dispersion 
of the gas in the cloud, the temperature of the CO-emitting 
gas, and the filling factor of bright CO emission ([Glover X 


Mac Low|2011|[Shetty et al.|2011a|b| [Narayanan X Hopkins 


2013). Increasing the gas temperature increases the typical 
brightness temperature of the emitting gas, thereby increas¬ 
ing Wco and decreasing Xco- However, as a change in the 
temperature also affects the excitation of the CO molecules, 
the scaling of Wco with tempera ture is not linear, bu t in¬ 
stead is closer to Wco cx: (Shetty et al. 2011b). In¬ 


creasing the velocity dispersion also increases Wco by in¬ 
creasing the CO linewidth, but the increased velocity dis¬ 
persion can also lead to lower CO line opacities and hence 
lower brightness temperatures, so again the dependence is 
sub-linear ( Shetty et al.|20lTb ). Finally, Wco is sensitive to 
the fraction of the cloud traced by bright CO emission, as 
in cases where this is small, the mean integrated intensity 
is much smaller than the peak integrated intensity, owing to 
the effects of beam dilution. In solar metallicity clouds situ¬ 
ated in a standard Galactic radiation field, the filling factor 


of the CO emission is of order unity (see e.g. Wolfire, Hol- 
lenbach X McKee 2010 or Section 3 below), and so small 


variations in its value from cloud to cloud have little effect 
on Xco- At lower metallicities, however, the CO filling fac¬ 
tor can become very small ( Glover X Mac Low|2011 Glover 
X Clark||2012b ) and hence Xco can increase dramatically. 

If we also allow the H2 column density to vary, then this 
adds an additional parameter to the problem. However, in 
real clouds, it is unlikely that all of these parameters vary 
independently. In particular, we expect the filling factor of 
CO-bright gas to be sensitive to the mean extinction of the 
cloud ( Glover X Mac Low|2Qll ), which at fixed metallicity is 
directly related to Xh 2 - Nevertheless, in order to understand 
why Xco is close to constant in local GMCs, we need to 
understand why these four parameters do not all vary by 
large amounts. 

If GMCs are in virial equilibrium with linewidths that 
satisfy the standard Galactic size-linewidth relationship 


(Solomon et al.|1987||Scoville et al.|1987||Roman-Duval et al. 

20T0I 


D \ 

a ~ 0.7 ( -- ) kms“^, 

IpcJ 


( 2 ) 


where a is the linewidth and R is a measure of the char¬ 
acteristic size of the cloud, then we would expect that Xh 2 
would be approximately constant from cloud to cloud (Lar¬ 


son 1981) and that the CO linewidth would scale only weakly 


with the cloud mass as a oc (Bolatto, Wolfire X Leroy 


2013), thereby explaining much of the constancy in Xco- 


However, it has become clear over the last few years that 
this is not the only possible explanation and that virial equi¬ 
librium is not required in order to explain a constant Xco- 
The reason is that even if GMCs are not virialized, their 
mean H2 column densities are unlikely to vary by a large 
amount. Clouds with low values of Xh 2 will provide insuffi¬ 
cient dust shielding to allow widespread CO formation and 


hence will not be identified as “molecular” clouds (Glover X 
Mac Low 2011 Clark X Glover|2014 ), while clouds with high 
Xh 2 will form stars rapidly, and the resulting stellar feed¬ 
back will prevent their column densities from growing too 


large ( Feldmann et al.||2012 Narayanan X Hopkins||2013 ). 

The small variation that we find in the mean CO bright¬ 
ness temperature of local GMCs can also be easily under¬ 
stood. Since ^^CO emission from GMCs is usually optically 
thick, the observed brightness temperature is determined by 
the density and temperature of the gas at the point in the 
cloud where r — 1 ( Bolatto, Wolfire X Lero;^|2013 ). These 
values vary much less from cloud to cloud than do the mean 


values for the cloud as a whole ( 

Molina, Glover X Feder- 

rath||2011 Narayanan X Hopkins 

2013), and so there is lit- 


tie variation in the resulting mean brightness temperatures 
( Roman-Duval et al.||2010 Narayanan X Hopkins 1 2 OTS ). 


Putting all of this together, we see that we can under¬ 
stand fairly well why Xco is approximately constant for 
local GMCs. However, the behaviour of Xco as we move to 
environments with much stronger UV radiation fields and 
higher cosmic ray fluxes is less easy to predict. We would 
expect that as we increase the strength of the ambient ra¬ 
diation field, more photodissociation of CO will occur. This 
will decrease the filling factor of the CO emission and may 
also decrease the CO linewidth if the CO-emitting gas is no 
longer well-distributed throughout the volume of the cloud. 
These effects will tend to decrease Wco and hence increase 
Xco- At the same time, we also expect a stronger UV field 
to lead to greater dissociation of H2, and hence a decrease in 
Xh 2 - In addition, we would expect the cloud to be system¬ 
atically warmer, owing to the higher photoelectric heating 
rate in low Ay regions, and the greater cosmic ray heating 
rate in high Ay regions. This means that the mean bright¬ 
ness temperature of the CO emission may also be larger. 
Together, these two effects will tend to decrease Xco- It is 
not immediately clear which set of effects will dominate, and 
hence whether one would expect Xco to increase or decrease 
as the strength of the UV radiation field and the cosmic ray 
flux increase. 


In this paper, we investigate this issue with the help of 
numerical simulations. We simulate the coupled chemical, 
thermal and dynamical evolution of several model clouds in 
a variety of different environments, spanning a range of 100 
in terms of UV field strength and cosmic ray ionisation rate. 
We model and analyze the resulting emission in the 1-0 line 
of ^^CO, and thereby are able to explore how Xco varies as 
a function of the external environment in clouds with a fixed 
mass and mean volume density. In Section 2, we briefly de¬ 
scribe our numerical approach, and in Section 3, we explore 
how the changing environment affects the physical structure 
of the clouds. In Section 4, we examine how Xco varies from 
cloud to cloud, and how this depends on the definition used 


© 0000 RAS, MNRAS 000 , 000-000 



























































Does Xqo depend on the SFR ? 3 




total column density [cm Hg column density [cm 


0.1 


1.0 

Wco [K 


10.0 

km s“‘ 


100.0 


Figure 1. The columns show the total column density (left), the H 2 column density (middle) and the velocity-integrated CO intensity 
in the 1-0 line (right) for the “standard” SFRq (top row), 10X SFRq (middle row) and 100x SFRq (bottom row) in the 10^ M© clouds. 
All of the images were made at the point at which the density exceeded 10® cm“^ in the first collapsing core in each simulation. 
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Figure 2. As Figurebut for the 10^ Mq clouds. Note the change of physical spatial scale, as illustrated in the top left panel. 


for the cloud boundary. We discuss some implications of our 
results in Section 5 and conclude in Section 6. 


2 NUMERICAL METHOD 
2.1 Details of the algorithms 

The computations presented in this paper were performed 
using a modified version of the publicly available smoothed 


particle hydrodynamics (SPH) code GADGET-2 (Springel 
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2005). A major modification to the original GADGET-2 is 


the inclusion of time-dependent chemistry. The version we 
use here follows the formation and destruction of H2 as in¬ 
troduced by Glover & Mac Low ( 2007a|b ) in addition to the 
simplified treatment of the GO chemistry that was proposed 
by Nelson & Langer (1999). In Glover & Glark (2012a), it 


was shown that this GO network has a similar accuracy to 


the more exhaustive treatment of Glover et al. (2010), but 


incurs around one third of the computational cost. We do 
not include freeze-out of GO onto dust grains in this model, 
but as this generally occurs only in regions where the ^^GO 


emission is already optically thick ( Goldsmith||2001 ), we do 
not expect this omission to significantly affect our results. 

As well as following the chemical evolution of the gas, 
we also model its thermal evolution. We account for dynam¬ 
ical heating due to shocks and adiabatic compression and 
cooling due to adiabatic rarefaction in the same fashion as 
in the unmodified version of GADGET-2. In addition, we 
also account for the main radiative and chemical heating and 
cooling processes occurring in the ISM. These include fine 
structure cooling from G^, G and O, molecular line cool¬ 
ing from H2 and GO, photoelectric heating and cosmic ray 
heating. At high gas densities, collisions between gas parti¬ 
cles and dust grains also play an important role in regulating 
the thermal energy balance, cooling the gas if Tgas > Tdust 
and heating it if Tdust > Tgas- Eull details of how we treat 
these processes and a number of other, less important, con¬ 
tributors to the overall thermal energy balance can be found 
in our previous papers ( Glover Mac Low|20Q7gt|b Glover 


et al.|[MT0 Glover Glark||2012a ), and a summary of the 

most important processes included can be found in Eigure 
4 in [Glark et al. (2013). 

The attenuation of the interstellar radiation field 
(ISRE) is treated using the TreeGol algorithm, introduced 
by [Glark, Glover &: Kless^ (2012). In this paper, the spec¬ 
tral shape of the ISRE is based on the prescription of [Draine| 
(1978) in the ultraviolet and Black (1994) at longer wave¬ 


lengths. The strength of the ISRE is varied in the different 
simulations as described in Section l231 below. The clouds are 
assumed to be bathed in a uniform ISRE, and TreeGol is 
used to compute the attenuated spectrum that reaches each 
SPH particle in the computational volume. 

Since the above-mentioned papers, we have made two 
significant changes to the chemical model. The first was to 
update our treatment of the photodissociation of GO from 


the prescription given in Lee et al. (1996), to that described 


in the recent paper of Visser, van Dishoeck, & Black (2009). 
The second change that we have made to the chemical model 
is the inclusion of the effects of cosmic-ray ionisation of 
atomic carbon 

C + c.r.C++e“, (3) 

and cosmic-ray induced photodissociation of G and GO 
C + 7cr G"^ -h e , (4) 

GO-h7cr ^ G-hO. (5) 

These processes were not included in the chemical model 


described in Glover & Glark (2012a). We assume that the 


rates of all three processes are proportional to Ch , the cosmic 
ray ionisation rate of atomic hydrogen. Eor processes and 
1^ we use the rates given in [Woodall et al. (2007) as a basis, 
but rescale them by a factor Ch/Ch,wo7, where Ch,wo7 is the 


value of the cosmic ray ionisation rate adopted by [Woodall[ 
et al. (2007). Similarly, to compute the rate of process ^ we 


use the value given in Gredel, Lepp & Dalgarno (1987) as a 


basis, but rescale it to make it consistent with our choice of 
Ch. 

The ^^GO (1-0) emission maps that form the basis of the 
analysis in this paper were created using the RADMG-3D 
radiative transfer cod^ The level populations in the non- 
LTE limit are computed using the the large velocity gradient 
approximation ( Sobolev|1957 ), as implemented in RADMG- 
3D by Shetty et al. ^ 2011a|b ). We assume that collisions with 
H2 dominate the excitation of the GO rotational levels and 
adopt the values for the collisional excitation rate coefficients 


given in the Leiden Atomic and Molecular database (Schoier 
et al.[[2005 ). 


2.2 Cloud properties 

All of the clouds in this study start as uniform spheres of 
gas, onto which a three-dimensional turbulent velocity field 
has been superimposed. Two cloud masses are examined in 
this paper: a “low-mass” cloud, with a mass of 10^ Mq, and 
a “high-mass” cloud of mass 10^ M©. In most of the runs, 
we take the initial hydrogen nuclei number density of the 
gas to be no = 100 cm“^, yielding initial radii for the low- 
and high-mass clouds of 8.8 pc and 19 pc, respectively. In 
addition, we performed two simulations of high-mass clouds 
with a much higher initial density, no = 10^ cm“^; these 
denser clouds had an initial radius of around 4 pc. 

The mass resolution in this study is kept fixed, with the 
mass of an SPH particle being set at 0.005 M©. The low- 
mass clouds therefore have 2 x 10® SPH particles, while the 
high-mass clouds have 2 x 10^ SPH particles. The minimum 
resolvable self-gravitating mass-element in this calculation 
is therefore 0.5M© in both cases ( Bate Burkert|1997 ). In 
practice, this means that we can follow the collapse of the 
cloud until the gas number density reaches a value of around 
10 ®cm“^. At this point, we halt the collapse and perform 
the analysis that is presented in this paper. We note that we 
would not expect the results presented here to change if we 
were to follow the collapse to higher densities, as at these 
densities, the GO (1-0) line will be highly optically thick (see 
e.g. Goldsmith|[2001[ ), and in any case much of the GO will 
be frozen out onto dust grains. 

The velocity fields are generated with a ‘natural’ mix of 
solenoidal and compressive modes (i.e. a ratio of 2:1). This 
is generated on a 128^ grid. This velocity field is left to freely 
decay in shocks, rather than being continuously driven dur¬ 
ing the course of the simulation. In most of the simulations, 
the initial energy in the velocity field, Ekin, was set to be 
half the initial gravitational potential energy of the cloud, 
Egrav, so that the clouds are initially in virial equilibrium. 
This means that the initial RMS turbulent velocity I’rms is 
2.4 kms“^ for the low-mass clouds and 5.2 kms“^ for the 
high-mass clouds. In terms of the virial parameter 


Ekii 

Egra 


( 6 ) 


^ http://www.ita.uni-heidelberg.de/~dullemond/software/radmc- 
3d/ 
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these clouds have a^iv = 0.5. In addition, we also examine 
one high density cloud for which we set avir = 2, so that the 


cloud is initially unbound (see Section 4.3). Note that most 
of the simulations presented here adopt the same random 
seed for the turbulent velocity field. However we also run 
one of the simulations with a different seed, to gauge the 
sensitivity of our results to the underlying cloud structure 
that is created by the turbulent motions. 

From previous modelling, it has been found that in cold 
gas with densities of around 100 cm“^, most of the hydrogen 
is in the form of H 2 . In contrast, the carbon in this gas 
is still predominantly in the form of (see e.g. Clark et 
al. 2012 Smith et al. 2014). The initial chemical state of 


the gas in the majority of the runs presented in this study 
is motivated by this previous work. We start most of our 
runs with all of the hydrogen already in the form of H 2 , but 
assume that the carbon and oxygen are present in the form of 
and O, respectively. We adopt total abundances (relative 
to the number density of hydrogen nuclei) for the carbon 
and oxygen nuclei of = 1-4 x 10“^ and xo = 3.2 x 10“^ 
respectively, consistent with the values measured in the local 
ISM ( Sembach et al.|[2000 ). Similarly, we adopt a dust-to¬ 
gas mass ratio of 0.01, consistent with the value in the local 
ISM. 

In our runs with very high UV fields, it is unclear 
whether starting with all of the hydrogen in molecular form 
is a good approximation, as in this case, one would ex¬ 
pect the equilibrium H 2 fraction in low density gas to be 


much smaller than in the models of Clark et al. (2012) and 


Smith et al. (2014). To address this uncertainty, we ran two 


additional high-mass models with the hydrogen initially in 
atomic form. 

The post-processing of the simulation data in RADMC- 
3D first requires that the SPH particle data is interpolated 
onto a regular Cartesian grid. This is done using the stan¬ 
dard SPH smoothing formalism. To ensure that we catch 
small, high density pockets of gas, we employ a RADMC- 
3D cell-size of 0.068 pc, such that the high-mass cloud calcu¬ 
lations have 656^ cells and the low-mass cloud calculations 
have 256^ cells. We find that such a resolution is sufficient 
to get converged values for the probability density function 
of CO (1-0) emission, from which it follows that the vari¬ 
ous mean values that we examine later are also converged. 
Once the position-posit ion-velocity maps are obtained, we 
then integrate the emission along the z-axis to create maps 
of the integrated intensity, VFco(i-o)- 


2.3 The “star formation rate” 

In this study, we vary two of the environmental conditions 
that can affect the chemical balance of clouds: the strength 
of the ISRF and the cosmic-ray ionisation rate. Both of these 
are thought to vary with the local star formation rate, and 
so we assume here that the strength of these processes can 
be used as a proxy for the star formation rate. 

As mentioned above, the ISRF used here is taken to 


have a shape described by a combination of the Black (1994) 


and Draine (1978) radiation fields. In one set of runs - those 


representing clouds in an environment with a star formation 
rate similar to that in the local ISM, which we denote as 
SFRo - we adopt the same normalization for the ISRF as in 


the papers of Draine (1978) and Black (1994). Our fiducial 


ISRF therefore has a strength Go = 1.7 in Habing (1968) 
units. 

For this study, we are mainly interested in how the ISRF 
heats the gas and affects its chemical state. Photons with en¬ 
ergies above 6 eV are responsible for the photoelectric heat¬ 
ing (the dominant heat source in low extinction regions of 
these clouds), and photons with energies above 11.2 eV and 
11.5 eV are responsible for dissociating H 2 and CO respec¬ 
tively. Since most of the photons in this part of the ISRF 
come from massive, young stars, it is reasonable to assume 
that, to a first approximation, the strength of the relevant 
portion of the ISRF scales linearly with the local star for¬ 
mation rate. At longer wavelengths, the ISRF is dominated 
by older stellar populations and this assumption is less well- 
founded. However, the strength of the ISRF at these wave¬ 
lengths has little effect on the temperature or chemistry of 
the gas, and so for simplicity, we assume that in regions with 
higher star formation rates, we can simply scale the entire 
ISRF upwards, rather than changing its spectral shape. 

We also assume that the cosmic ray ionisation rate 
scales linearly with the star formation rate. This assump¬ 
tion is reasonable: supernova remnants are the main source 
of cosmic rays in the Galactic ISM, and the lifetime of a typ¬ 


ical cosmic ray within the Galaxy is around 15 Myr (Ferriere 


2001 ), so the cosmic ray energy density, and hence the cos¬ 


mic ray ionisation rate, should track the star formation rate 
fairly closely. 

In this study, in addition to the runs representing the 
behaviour of clouds in the local ISM, with SFR = SFRo, we 
also perform simulations where we increase the strength of 
the ISRF and the cosmic ray ionisation rate by factors of 10 
or 100, corresponding to star formation rates of 10 SFRo or 
100 SFRo, respectively. 


3 GENERAL CLOUD STRUCTURE 

In the left-hand columns of Figuresandwe show column 
density images of the gas, taken at a point in the simulations 
just before the onset of star formation in our fiducial clouds 
- i.e. those seeded with our standard turbulent velocity field, 
and which started life with all their hydrogen in the form 
of H 2 . This particular stage in the cloud’s evolution is when 
we perform the analysis that forms the basis of this paper. 
For most of our simulated clouds, this corresponds to a time 
t ~ 1-2 Myr after the beginning of the simulation, as shown 
in Table [ 1 ] 

In both the 10^ M© and 10^ Mq clouds we see a sim¬ 
ilar change in behaviour of the strength of the ISRF and 
cosmic ray ionisation rate are increased to mimic progres¬ 
sively higher rates of ambient star formation activity: the 
filamentary structure imposed by the cloud’s turbulence and 
self-gravity is less pronounced as we move to higher SFRs. 
This is simply a result of gas temperatures rising due to the 
combined effects of the enhanced photoelectric emission and 
cosmic-ray ionisation rates. This systematic increase in the 
temperature as we increase the SFR is also illustrated in Fig¬ 
ure where we show temperature-density phase diagrams 
for four of our runs, colour-coded by the GO abundance. 

In addition, we see from Figures and that the over¬ 
all morphology of the clouds changes as the SFR increases. 
In this respect, there are two trends. The first is that the 
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Table 1. List of simulations. All number densities are given with respect to the number of hydrogen nuclei, since this is invariant of the 
chemical state of the gas. The definition of avir can be found in Equation while SFRq is defined in Section |2.3| The property tend 
denotes the point at which the simulation is stopped, and the cloud details are fed to the line radiative transfer code. This is also the 
time at which the first star forms in the simulation. 


Mass (Mq) 

no (cm ^) 

«vir 

SFR (SFRo) 

Notes 

^end (Myr) 

10 ^ 

100 

0.5 

I 


1.83 

10"^ 

100 

0.5 

10 


2.09 

10^ 

100 

0.5 

100 


I.9I 

10"^ 

100 

0.5 

100 

Different seed 

2.18 

10® 

100 

0.5 

I 


1.17 

10® 

100 

0.5 

10 


1.52 

10® 

100 

0.5 

100 


1.39 

10® 

100 

0.5 

I 

Atomic ICs 

1.31 

10® 

100 

0.5 

100 

Atomic ICs 

1.26 

10® 

10^ 

0.5 

100 

Galactic Centre style cloud 

O.I 

10® 

0 

1— 1 

2 

100 

Galactic Centre style cloud 

O.I 


dense, filamentary structure becomes confined to a progres¬ 
sively smaller region as the SFR increases. This makes sense, 
as the gas is more easily structured near the centre of the 
cloud where the extinction is higher and the heating effects 
of the photoelectric emission are less. As the SFR increases, 
the heating by photoelectric emission becomes progressively 
more important in shaping the gas, and so the regions struc¬ 
tured by the turbulence and gravity are those towards the 
centre. 


The second effect is that the structure of the cloud en¬ 
velope changes. The SPH particles have open (i.e. vacuum) 
boundary conditions, and so the particles are free to expand 
into the surrounding space. The stronger the ISRF, the more 
the surface of the cloud is heated, causing it to “boil off” 
into the void, creating the halo that we can see around the 
clouds. In the case of the 100 SFRo clouds (i.e. those illu¬ 
minated by radiation fields and cosmic ray fluxes that are 
100 times greater than the values in the local ISM; see Sec¬ 
tion 2.3), the ISRF is so strong that this halo becomes hot 


enough to push back on the cloud, creating the sharp bound¬ 
ary that we see in Figures^ andAlthough dramatic, this 
feature of the increasing SFR is actually of little interest to 
this current study, as both the H 2 and CO lie within the 
cloud boundary. 


The middle panels in Figures and show the H 2 col¬ 
umn densities images at the same point in the cloud evolu¬ 
tion. In general, we see that the H 2 column density tracks 
the total column density extremely well for the bulk of the 
cloud. Given that all the clouds in these two figures were 
initialised with all their hydrogen in molecular form, this is 
not surprising. 


The main difference between the left and middle pan¬ 
els in Figures and is that the H 2 in the outer envelope 
has been photodissociated by the ISRF, and so does not 
trace the very low density (both column and volume) that 
surrounds the main cloud structure. We also see that the 
column density at which H 2 starts to appear increases with 
increasing SFR, moving from around I0^° cm“^ to around 
10^^ cm“^ as the SFR increases from the solar neighbour¬ 
hood value to a factor of 100 greater. This behaviour is con¬ 
sistent with the predictions of detailed models for the struc¬ 


ture of photodissociation regions (see e.g. Krumholz et al. 
2008 Sternberg et al.|20T4 ). The fact that a large change in 


the strength of the ISRF leads to a relatively small change 
in the column density of gas required to effectively shield H 2 
is a consequence of the fact that dust shielding is an expo¬ 
nential function of the column density, which implies that 
in the dust-dominated region, the critical column density 
scales only logarithmically with the SFR. 


The integrated emission in the I-O line of ^^CO, Wco, 
is shown in the right-hand column in Figures and and 
shows the most striking variation with increasing SFR. For 
the SFRo clouds, we see that the surface filling factor of 
the bright CO emission is very similar to that of the H 2 
gas. Despite this, it is not a good tracer of the physical 
structure of the cloud: some features seen in the emission 
map are far more blobby than their H 2 counterparts and 
others are missing entirely. This behaviour is expected (see 
e.g. [Ballesteros-Paredes Mac Low||2002| [Goodman et ah] 


2QQ9| [Shetty et ar|20IIa|b| |Beaumont et ar| 2QI3) and is 

simply a result of the large optical depth of the I-O line - 
the features that we see are the r ~ I surfaces, which are a 
function of CO column density, velocity, and temperature, 
and so often bear little resemblance to the structure we see 
in the H 2 column density images. 


As we move to higher SFRs, we see that the surface 
filling factor of the CO emission progressively decreases, a 
feature we see in all of the cloud models run in this study. In 
the case of the 100 SFRo clouds, the effect is extreme, with 
emission falling below IKkms“^ for much of the cloud’s 
interior - the sensitivity limit in many large-scale surveys 
- and even below O.IKkms”^ along lines of sight with H 2 
column densities as high as a few 10^^ cm“^. 


However, in the cases with the higher SFR, we see that 
the peak CO emission becomes brighter towards regions of 
higher column density. This is a consequence of increased 
heating in the cloud, which results in a higher overall ex¬ 
citation temperature of the CO. The broader line-width in 
the gas also helps to the lower the optical depth, allowing 
more of the emission from the line to escape ( Shetty et al.| 
20IIa|b ). Due to the higher spatial resolution in the repro¬ 


duced image (the true resolution is the same in all images), 
this effect is perhaps most easily seen in the 10^ M© clouds in 
Figure Here we see the brightening of the clump towards 
the left of the Wco(i-o) image, which goes from having a 
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CO abundance 


Figure 3. Phase diagrams for four of our simulated clouds, il¬ 
lustrating the temperature structure of the gas at the end of the 
simulation. The diagrams are colour-coded using the mean CO 
abundance at each point in temperature-density phase space. 


peak brightness of around SOKkms ^ for SFRo, to nearly 
80Kkms“^ in the case of 100 SFRo 


4 X-FACTOR VARIATIONS 


4.1 Computing the X-factor 


In the previous section, we saw that as the SFR in the 
ambient cloud environment is increased, the CO emission 
traces progressively less of the H 2 distribution. However we 
have also seen that the peak emission in the cloud is higher 


in the case of high SFR. Bolatto, Wolfire & Leroy (2013) 


have suggested that these two effects may effectively cancel 
out, resulting in an X-factor that does not depend on the 
cloud’s environment. In this scenario, the greater emission 
that comes from the dense bright peaks in the high SFR case 
makes up for the fact that much of the cloud is not emitting 
in CO at all, since the molecule is easily photodissociated. 
Using our self-consistent models for Nhs and Wco, we can 
now explore this idea. 

In what follows, we compare how Xco varies as we al¬ 
ter the environmental conditions in our suite of turbulent 
clouds. We calculate Xco for each cloud via 


Xco = 


{NH 2 ) 

(M/co> 


(7) 


where the averages are taken over a specified set of the pixels 
in the images from Figures and 

In our evaluation of Equation we take two ap¬ 
proaches. First, we examine the values of Xco obtained 
when we average over all of the pixels in the images. In 
our second approach, we instead restrict the calculation to 
include only those pixels with CO I-O integrated intensities 
above some threshold value, fUco,cut- These two approaches 
mimic different types of observations. The first approach 
is more appropriate when considering unresolved molecular 
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Figure 4. Mean CO-to-H 2 conversion factor, Xco 5 in each of 
our clouds, as a function of our SFR proxy (see Section |2.3| 
for details). Red and blue circles denote the results from the 
no = 100 cm“^ clouds with fully molecular and fully atomic ini¬ 
tial conditions respectively. The orange circles correspond to a 
cloud with no = 100 cm“^ run using a different random seed for 
the turbulent velocity field. The dark and light green circles de¬ 
note the no = 10"^ cm“^ clouds with, respectively, avir = 0.5 and 
2. In the top two panels, we compute Xco using integrated inten¬ 
sities and H 2 column densities only for those lines-of-sight with 
an integrated intensity above the specified threshold, IFco,cut- lu 
the bottom panel, we show the results that we obtain if we do not 
apply an integrated intensity threshold. In each case, we compute 
Xco by dividing the mean H 2 column density in the considered 
area by the mean value of VFco from the same area. The horizon¬ 
tal dashed line shows the canonical local value of Xco 5 und the 
shaded area indicates the typical scatter in the value of Xco de¬ 
rived observationally for the Milky Way and other nearby spiral 
galaxies (see|Bolatto, Wolfire & Leroy|20I3|. 
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clouds that could inhabit the beam in an extragalactic CO 
observation. In this case, the CO emission may not trace 
all of the H 2 gas in the beam, i.e. there is a ‘dark’ molec¬ 
ular component surrounding the CO-bright portion of the 
cloud, as recently highlighted by |Wolfire, Hollenbach Mc-| 


Kee 


(2010). The second approach, in which we only exam¬ 


ine pixels with emission above ITco,cut, is more applicable 
to studies of nearby, resolved molecular clouds, where it is 
common to define the extent of a cloud in terms of its CO 
emission, a definition which inevitably leads to one ignor¬ 
ing those regions without detectable CO emission (see e.g. 
Pineda, Caselli Goodman||2008 Lee et al.|2Q14 ). 


The mean values of Xco that we derive from the simula¬ 
tions, and how they scale with the SFR, are shown in Figure 
1^ The top two panels show the results that we obtain for 
two different CO detection thresholds, VFco,cut = 3Kkms“^ 
and Wco,cut = 0.3 Kkms“^, while the bottom panel shows 
the results in the case where we do not apply a detection 
threshold. We also indicate the canonical local value of the 
X-factor, Xco,gal (the horizontal dashed line) and illustrate 
the typical scatter in the value of Xco derived for the disks 
of nearby spiral galaxies (the shaded region), as summarized 
in Bolatto, Wolfire X Leroy| ( |2013| ). 


4.2 Clouds with SFR = SFRo 

We begin our study of the effects of the local environment 
on Xco by looking at the behaviour of clouds located in 
an environment with a star formation rate similar to the lo¬ 
cal ISM, i.e. with SFR = SFRo (see Section 2.3). The first 


point to note is that when we consider clouds in this envi¬ 
ronment that have properties similar to those of local molec¬ 
ular clouds (i.e. virial parameters close to unity and mean 
initial densities around 100 cm“^), we find that Xco does 
not appear to be sensitive to mass. Both our low-mass and 
our high-mass clouds yield a mean X-factor of between 4-5 
X10^° cm“^K“^ km“^ s, consistent with the large body of 
literature suggesting that the X-factor is not strongly de¬ 
pendent on the cloud mass (see e.g. [Solomon et ^|1987| 
Blitz et al.||2007l [Hughes et al.||2010| [Bolatto, Wolfire X| 
Leroy|20T3 ). Also, we hnd that the X-factor is relatively in¬ 
sensitive to the value of Wco,cut, indicating that the CO 
emission is a good proxy for the total H 2 present in such 
clouds, even though the emission is clearly not tracing all of 
the molecular gas. 

While this is encouraging, it should also be noted that 
the values we recover for Xco in this case are roughly twice 
as large as the canonical Galactic value, Xco,gai, and on the 
edge of the range of values found for the disks of nearby 
spiral galaxies ( Bolatto, Wolfire X Leroy[[2013 ). The reason 
for this discrepancy is not completely clear, but there are 
several factors that may contribute to our recovering a sys¬ 
tematically higher Xco- The first of these stems from the 
way in which we treat the effects of self-shielding in our 
simulations. When calculating each particle’s sky maps of 
H 2 and CO column densities, the TreeCol algorithm in¬ 
cludes contributions from all the molecules present in the 
cloud, regardless of the velocity at which they are moving 
with respect to the target particle. In reality, much of the 
cloud will not contribute to the self-shielding of a given lo¬ 
cation, since the cloud’s supersonic motions will Doppler 
shift the molecules out of the line-profile. Neglecting this ef¬ 


fect means that we overestimate the effectiveness of H 2 and 
CO self-shielding, and hence overproduce the abundances of 
both molecules. We expect the error to be larger for H 2 than 
for CO, since self-shielding plays a much greater role in de¬ 
termining the H 2 abundance than the CO abundance, and 
therefore it will lead to a systematic increase in our values of 
Xco compared to those for real clouds. The importance of 
this effect will depend on the size of the velocity dispersion 
in the cloud, but should not depend strongly on our choice 
of SFR, and it therefore represents a systematic error that 
biases all of our derived values of Xco to slightly higher 
values. 

Another factor that may contribute to the discrepancy 
between our results and Xco,gai is our choice to start the 
simulations with the hydrogen in fully molecular form. Al¬ 
though this approximates the behaviour that we find in pre¬ 


vious studies of cloud formation (Clark et al. 2012 Smith 


et al. 2014), in practice the hydrogen in the assembling 
clouds is never completely molecular - there is always some 
atomic component associated with the cloud. In equilibrium. 


we expect the surface density of atomic hyd 
with the clouds to be of order 10 M© pc“^ ( 

rogen associated 

[Krumholz et al. 

2008 

[Wolfire, Hollenbach & McKee[[2010| | 

pternberg et al. 

2014 

), but the time taken to reach this equilibrium can be 


longer than the lifetime of the cloud. In our simulations with 
fully molecular initial conditions, we will therefore overesti¬ 
mate the H 2 content of the clouds by around 10-20%, and 
hence will overestimate Xco by a similar factor. In order 
to verify that the effect on Xco is not greater than this, 
we ran two simulations of our high-mass cloud in which 
we started with all of the hydrogen in atomic form. In one 
of these simulations, we adopted the solar neighbourhood 
value for SFR (i.e. SFR = SFRo), while in the other we set 
SFR = 100SFRo. The values of Xco that we derived from 
these two runs are indicated in Figure by the blue cir¬ 
cles. We see that, as expected, the use of atomic rather than 
molecular initial conditions results in a decrease in Xco, but 
that the effect is not large and seems to be independent of 
our choice of SFR. 

Finally, our decision to focus on clouds with a virial 
parameter avir = 0.5 may also lead to us deriving a larger 
value of Xco than that found for large samples of Galactic 
GMCs. As discussed in the introduction, we know that the 
mean integrated CO intensity of a CMC depends sensitively 


on the CO line-width for that CMC (see e.g. Shetty et al. 
2011a|b |. Increasing the line-width while keeping all other 
cloud properties the same leads to an increase in Wco and 
hence a decrease in Xco • Similarly, decreasing the line-width 
decreases Wco and increases Xco- This implies that Xco 
should depend to some extent on avir, in the sense that 
clouds with larger virial parameters that survive for long 
enough to form CO should have smaller values of Xco than 
similar clouds with smaller values of avir- In any large sam¬ 
ple of GMCs, we would not expect to find many clouds with 
Qvir 1: any clouds that do initially have F^kin F^grav 
will undergo gravitational collapse, inducing motions that 
rapidly increase F^kin driving a^iv towards 1 if the clouds 
continue to collapse or 0.5 if they settle into virial equilib¬ 
rium ( Ballesteros-Paredes et al.[[2011 ). On the other hand. 


there is both theoretical (Dobbs, Burkert X Pringle 2011 
Ward, Wadsley X Sills|20T4 ) and observational (Heyer et al 

2009) evidence that many GMCs are gravitationally un- 
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bound, with a^ir > 1. It is therefore plausible that the mean 
value of Qvir for the population of GMCs in the Milky Way is 
of order unity or larger. If so, then it follows that the clouds 
that we model in this study have systematically smaller CO 
line-widths than the average Galactic CMC, helping to ex¬ 
plain why the values of Xco that we derive for them are 
larger than Xco,gai. 


4.3 Dependence on the SFR 


We now come to the question at the heart of this paper: 
does the value of Xco vary with the local ambient SFR? 
While the results displayed in Figure clearly show that 
Xco is indeed a function of SFR, the extent to which the 
X-factor varies depends on several parameters, and we will 
now discuss these in turn. 

First, examining the red circles in Figure]^ we see that, 
in general, the mean X-factor derived from the cloud maps 
increases as SFR increases. Given that CO is more easily 
dissociated than H 2 , this makes sense: the CO is doing a 
progressively poorer job of tracing the true extent of the 
molecular gas as the ambient SFR increases. However we 
see that the effect is more pronounced for the high-mass 
clouds than for the low-mass clouds. Indeed, there is a slight 
decrease in the X-factor of the 100 SFRo case with respect 
to the 10 SFRo case for the low-mass clouds. So while the X- 
factor at solar neighbourhood values of the SFR (i.e. SFRo) 
is independent of mass, this is not true when the SFR is 
increased. 

The reason behind this behaviour in the low-mass 
clouds can be inferred from Figure where we show the 
fraction of the cloud’s total CO emission that comes from 
pixels below a given Wco- For progressively higher SFR, we 
see that fraction of the overall integrated emission coming 
from very bright lines of sight depend on the cloud mass. 
For the low-mass clouds at high SFR, the majority of the 
emission comes from regions with integrated intensities of 
10 Kkms“^ or greater, while the opposite is true for the 
high-mass clouds, in that they gain most of their emission 
from lines of sight with lower Wco- The idea summarised 


by Bolatto, Wolfire X Leroy (2013) that a few bright lines 


of sight can compensate for the lack of CO elsewhere in the 
cloud only appears to hold for our low-mass clouds. 

The reason why the clouds behave differently has to 
do with their column density distributions. In the low-mass 
clouds, the mean column density is low enough that CO is ef¬ 
ficiently photodissociated throughout much of their volume 
as the SFR is increased. The CO that survives is confined 
to the highly overdense regions produced by turbulent com¬ 
pression and gravitational collapse, in which the shielding is 
more effective. Because of the high density, the excitation 
temperature of the CO in these regions is high, and so is 
the integrated intensity of the CO emission. These clouds 
therefore have many regions with very low Wco and a few 
with very high Wco that dominate the total emission. In 
the high-mass clouds, on the other hand, the mean column 
density of the gas is larger and the CO is therefore better 
able to resist photodissociation. Because of this, CO survives 
in lower density regions in these clouds even when the SFR 
is increased significantly. The excitation temperature of the 
CO in these regions is lower than in the dense cores, and so 
the integrated CO intensities coming from the lower density 
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Figure 5. Fraction of the total emission that arises from pixels 
of a given Wco and below. 


gas are smaller, but the much larger filling factor of these 
regions compensates for this, and leads to their contribution 
dominating the total emission from the cloud. 

It is plausible that molecular clouds in environments 
with higher SFRs will be denser, on account of the 
higher ambient gas pressure in these environments (see e.g. 


Elmegreen, Klessen X Wilson 2008 Colombo et al. 2014). 


Some support for this picture comes from observations of 
molecular clouds in the Central Molecular Zone of the Milky 
Way. The star formation rate in this region is roughly 50-100 
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times that in the solar neighbourhood (Bonatto & Bica 


2011 


Longmore et ah 2013a) and the molecular clouds located 
there have high mean densities, of the order of 10^cm“^ 


more (see e.g. Dahmen et ah 1998[ Longmore et ah 2012 


2013b| ). There is also evidence to suggest that these clouds 


have higher values of avir than local clouds (|Kauffmann, 
Pillai Zhang||2013 Johnston et al.||2014 ). 


To investigate this scenario, we performed two addi¬ 
tional simulations that probe this regime. We modelled 
two clouds with mass 10^ Mq and initial number density 
no = 10^ cm“^, in which the hydrogen was taken to be 
initially in molecular form. One of the clouds was initial¬ 
ized with a turbulent field resulting in a virial parameter 
Q^vir = 0.5 (i.e. virial equilibrium), while the other had an 
Qvir = 2 (i.e. the cloud was gravitationally unbound). These 
are shown as the green circles in Figure We see that the X- 
factors for these clouds are now much closer to the standard 
Milky Way value, demonstrating that the physical proper¬ 
ties of the clouds play a significant role in controlling the 
X-factor. We also see that, as expected, the X-factor is lower 
in the case where the amount of turbulent kinetic energy in 
the cloud is higher. 


Another factor that strongly affects the mean X-factor 
derived from the cloud maps is the value used for deciding 
whether or not a pixel should be included in the averaging 
process, Wco,cut- If we set Wco,cut = 0, and hence include 
all of the lines of sight in our synthetic emission maps when 
determining (iVH 2 ) and (Wco), then we find that Xco in¬ 
creases by between a factor of three (for the low-mass clouds) 
and an order of magnitude (for the high-mass clouds) as we 
increase the SFR from SFRo to 100 SFRq. However, if we 
restrict our calculation of the X-factor to only those lines of 
sight with detectable CO emission, we find that the depen¬ 
dence of Xco on the SFR becomes much weaker. For exam¬ 
ple, if we set Wco,cut = 3 K km s“^, we find that the X-factor 
increases at most by a factor of two as we increase the SFR 
from SFRo to 100 SFRq. This suggests that the main effect 
driving the increase in Xco with SFR in the Wco,cut = 0 
case is a large increase in the number of lines of sight that 
have molecular hydrogen but no detectable CO emission, 
or in other words that pass through “CO-dark” molecular 
gas ( Wolfire, Hollenbach X McKee|2010 ). In those regions of 
the clouds that remain CO-bright, the relationship between 
CO emission and H 2 column density varies only weakly with 
the SFR. This behaviour is similar to what is seen in obser¬ 
vational studies of the dependence of Xco on metallicity: 
observations that focus only on CO-bright regions find a 


weak dependence of Xco on Z (e.g. Wilson 1995 Bolatto 


et al. 2008), while those that sample both CO-bright and 


CO-dark gas find a much stronger dependence (e.g. Israel 


1997 Bolatto et al. 2011). Finally, it is also interesting to 


note that the X-factor for the high density, 100 SFRo clouds 
is insensitive to the value of Wco, cut- This is because these 
clouds are so dense that the “skin” in Ay-space that contains 
H 2 but not CO is so thin that very little mass is contained 
within it. As such, the CO typically always does a good job 
of tracing the molecular content of these clouds. 


5 DISCUSSION 


5.1 Consequences for the Kennicutt-Schmidt 
relation 

The results presented in the previous section demonstrate 
that if we restrict our attention to those portions of molecu¬ 
lar clouds that are traced by relatively bright CO emission, 
then the X-factor that we have to adopt in order to correctly 
convert from the CO intensity of these regions to their H 2 
column density does not vary strongly with the local envi¬ 
ronment. Increasing the size of the assumed star formation 
rate (and hence the strength of the ISRF and the cosmic 
ray ionisation rate) by a factor of 100 leads to an increase in 
Aco of around a factor of two at most. This result is con¬ 


sistent with the finding by Pineda et al. (2009) and Hughes 


et al. (2010) that in the LMC Aco does not vary strongly 
with the strength of the ISRF, since these authors deter¬ 
mine the masses of the clouds that they study using virial 
mass estimates that are sensitive only to the properties of 
the CO-bright gas. 

Unfortunately, the number of galaxies for which we 
can distinguish between CO-bright and CO-faint portions 
of GMCs is very small. This is routinely done in the Milky 
Way, has already been done to some extent in the Magellanic 
Clouds, and with the advent of ALMA becomes possible to 
do in more distant members of the Local Group. However, in 
most extragalactic observations, the best spatial resolution 
that we can achieve is comparable to or larger than the size 
of the individual clouds (see e.g. Colombo et al. 2014). In 
this case, we cannot easily distinguish between CO-bright 
and CO-faint portions of the clouds, and so the results that 
we present for the case when Wco, cut = 0 give the best guide 
to the behaviour of Aco in these systems. 

For these unresolved clouds, our results show that Aco 
increases substantially as we increase the SFR, provided that 
the other properties of the clouds remain unchanged (i.e. 
provided that their characteristic densities or velocity dis¬ 
persions do not also vary as functions of the SFR). One im¬ 
portant implication of this is that in order to be able to use 
CO emission as a reliable tracer of molecular mass, we need 
to understand the nature of the underlying cloud population 
and how this varies as a function of the star formation rate. 

If we make the simplifying assumption that cloud prop¬ 
erties do not depend to any great degree on the local value 
of the SFR, then our results imply that the standard in¬ 
terpretation of the Kennicutt-Schmidt relation may need to 
be reassessed. If the observed CO luminosity is dominated 
by high-mass clouds that have number densities similar to 
those of our fiducial cloud models, then our simulations sug¬ 
gest that Aco should scale with the ambient rate of star 
formation as roughl}!^ 

Xco oc (8) 

Provided that the SFR does not vary widely over the area 
of the observational beam, it then follows that 

Aco oc (9) 

Measurements of the Kennicutt-Schmidt relation at gas sur- 


^ Note that this scaling is approximate: it is not a precise fit to 
the simulation results. 
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face densities where the molecular component dominates 


e.g. Bigiel et'8dT]|2008 Blanc et'ar]|20Q9 Bigiel et'ar]|2011 


|Liu et al.||2011 ) find a power-law relationship that we can 

write as 


SSFR « 


-^obs 


( 10 ) 


where A^obs denotes the observed power-law index. In prac¬ 
tice, however, Emoi is usually not measured directly; instead, 
it is inferred from the CO luminosity through the applica¬ 
tion of some constant CO-to-H 2 conversion factor, which in 
normal spiral galaxies is generally taken to have the canoni¬ 
cal Galactic value. Therefore, we can more accurately write 
the above relationship as 


^SFK OC (i;cO,em-^CO,gal)^°^® , 


( 11 ) 


where Sco,em is the surface density of CO emission. Now, 
if the actual CO-to-H 2 conversion factor is not fixed, but 
instead varies with Esfr according to Equation then it 
follows that the actual relationship between Esfr and Emoi 
can be written as 


SsFR (X 


A^obs/2 

SFR 


( 12 ) 


Rearranging this, we find that after correcting for the de¬ 
pendence of Xco on the star formation rate, the true rela¬ 
tionship between Esfr and Smoi can be written as 


SsFR « sir 


where 


-Nact — 


2 iVobs 

2 T Nobs 


(13) 


(14) 


Recent work by Bigiel et al. (2008) has suggested a value 


of Nobs = 1 from a sample of nearby galaxies. In this case, 
the True’ index of the Kennicutt-Schmidt relation would 
be around Nact = 2/3, implying that CO is a progressively 
worse tracer of the star formation rate as one moves to more 


extreme environments (see the discussion in Shetty et al. 
2014b| ). However, more recently, there have also been claims 
that Nobs < 1 (i-c. that the Kennicutt-Schmidt relation is 
sub-linear), even when assuming a constant X-factor (e.g. 
Ford et ^|2013| [Shetty et al.||2Q13 2014a), with the most 
likely value lying at around 0.76. In this case, we would 
predict a value of Nact = 0.55. 

It is important to note that this conclusion depends 
strongly on our assumption that the cloud properties do not 
vary strongly as a function of the local star formation rate. 
If this is not the case, and clouds in regions with higher 
values of Ssfr are systematically denser and/or more tur¬ 
bulent than clouds in less actively star-forming regions, then 
it is plausible that the decrease in Xco caused by the higher 
densities and turbulent velocities could be large enough to 
overwhelm the increase caused by the higher star formation 
rate. Indeed, there is evidence that this is the case in the 
centres of many spiral galaxies: measurements of Xco often 
find that it decreases close to the centre (see e.g. Sandstrom 
eTZIpoBl l, while we know from the study of the Central 
Molecular Zone of our own galaxy that molecular clouds in 
this environment are very dense and highly turbulent (see 
e.g. the summary in the review of [Molinari et ah |2014| ). 
Similarly, observational determinations of Xco in nearby 
ultraluminous infra-red galaxies (ULIRGS) or rapidly star¬ 
forming systems at high redshift typically find values that 


are smaller than the canonical local value (see e.g. [Downes X 
' Solomon||1998| [M agdis et al.|2011[ [Ostriker X Slietty||2Qll| 


Hodge et al. |2012| |Fu et al.| 2013| or the recent review by 
Casey, Narayanan Cooray|2014| . However, these observa¬ 
tions probe regions with much higher gas surface densities 
and turbulent velocities than are found in the local ISM, 
once again suggesting that the decrease in Xco caused by 
the increased densities and velocities is able to overcome 
the increase caused by the higher radiation field strength 
and cosmic ray ionization rate. 

Ultimately, therefore, the real message to take away 
from this study is that it is necessary to understand the 
nature of the cloud population before one can properly in¬ 
terpret the CO emission. Unless one understands how the 
properties of the clouds depend on the local star formation 
rate, it is impossible to be certain whether Xco increases 
with increasing Ssfr (as is the case in our simulations at 
fixed ttvir), decreases with increasing Esfr (as appears to 
be the case in the centres of many spirals), or remains ap¬ 
proximately constant. Further compounding this problem is 
the fact that CO measurements are one of the main observa¬ 
tional tools used for determining the mass of molecular gas 
in clouds, potentially leaving us with a circular argument. 
Supplementary information on cloud masses, from tracers 
such as dust emission, can probably help to break this cir¬ 
cularity, but a detailed examination of how best to do this 
is outside of the scope of this paper. 

5.2 Comparison with previous work 

It is interesting to compare our results on the dependence 
of the X-factor on the SFR with those of previous numerical 


studies that have looked at this problem. Bell et al. (2006) 


use the time-dependent UCL_PDR code (Papadopoulos et al. 
2002 Bell et al.|2005 ) to examine how Nco varies as a func¬ 


tion of depth within a semi-infinite slab of gas as they vary 
a number of different physical parameters, including the ra¬ 
diation field strength. They find that at high ^v, increasing 
the radiation field strength has little effect on Nco: changing 
Go by a factor of 10® leads to no more than a factor of three 
change in Nco- At low ^v, on the other hand, the depen¬ 
dence of Nco on the radiation field is much stronger. At very 
low Hv, the gas is CO-faint and Nco is very large, but as we 
move to higher Ay, the CO content of the gas grows, and the 
value of Nco approaches the canonical Galactic value. The 
depth into the cloud at which this transition occurs depends 
on the radiation field strength, with stronger fields implying 
that the transition occurs at higher Ay . It is not straightfor¬ 
ward to convert from these results to a cloud-averaged Nco, 
but it is clear that they are at least qualitatively consistent 
with our findings that Nco measured only for CO-bright 
regions does not vary strongly with the SFR, while Nco av¬ 
eraged over the whole cloud (including the CO-faint regions) 
shows a much stronger dependence. 


Feldmann et al. (2012) have also explored the depen¬ 


dence of Nco on the local radiation field. They performed a 
large hydrodynamical simulation of the formation of a typ¬ 
ical L* galaxy (see also Gnedin N Kravtsov|2Qll ) and then 
post-processed the results of this simulation using a sub¬ 
grid prescription for Nco that assumes that it is primarily 
determined by the mean visual extinction of the individual 
clouds and the strength of the ISRF. For Go = 1, |Fdd^ 
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mann et al. (2012) use the results of Glover & Mac Low 
( |2011| to calibrate their subgrid model, while for Go ^ 1 
and Go ^ 1, they estimate the CO abundance (and hence 
Xco) using a simple remapping procedure that assumes that 
the clouds are in photodissociation equilibrium. They find 
that for clouds with mean extinctions Ay ~ 6 and below, 
the value of Xco depends strongly on the radiation field 
strength. For higher column density clouds, on the other 
hand, Xco becomes largely independent of the radiation 
field strength. As most of the molecular clouds that form 
in their simulations have relatively large column densities, 
they find that when you average over the whole population 
of clouds, the dependence of Xco on the strength of the 
ISRF becomes very weak. 

At first sight, these results would appear to contradict 
our finding that Aco can vary significantly as we change the 
SFR. However, there are a couple of important points that 
one should bear in mind. First, Feldmann et al. (2012) do 


not vary the cosmic ray ionisation rate, only the radiation 
field strength. Therefore, their models do not account for the 
increased destruction of CO by dissociative charge transfer 
with He^ that occurs when the cosmic ray ionisation rate 
is large. At high Ay, this is the main CO destruction mech¬ 
anism, and its inclusion in our models is one of the main 
reasons why we see at least some dependence of Aco on 
SFR even when we restrict our attention to portions of the 
cloud that are highly shielded and CO-bright. Second, the 
hydrodynamical simulations used as a basis for the |Feld-| 
mann et al. (2012) analysis have a resolution of only 60 pc 
and hence are sensitive only to the largest, densest clouds, 
which are naturally the ones least affected by changes in the 
SFR. 

Another large-scale numerical study of the dependence 
of Aco on environment was presented by [Narayanan et al.| 
(2012). They post-processed an extensive series of SPH sim¬ 
ulations of isolated and merging galaxies in order to de¬ 
termine the local chemical and thermal state of the gas 


and the consequent CO emission (see also Krumholz, Leroy 


A McKee 2011 and Narayanan et al. 2011). They found 


that in their models, Aco scales inversely with the molec¬ 
ular gas surface density as Aco oc In these models, 

Dsfr oc by construction, and so the implication is that 
Aco should decrease weakly with increasing star formation 
rate, scaling as Aco oc However, as with many of the 

observational studies of Aco in extreme environments men¬ 
tioned in the previous section, direct comparison of these 
results with our own is complicated by the fact that the 
mean properties of the clouds in regions with high SFR in 
these simulations are not the same as those of the clouds in 
regions with low SFR. 


Lagos et al. (2012) also studied the effects of changes 


in the strength of the ISRF and the cosmic ray ionization 
rate on the value of Aco using a simple semi-analytical ap¬ 
proach. They used the GALFORM semi-analytical galaxy for¬ 
mation model ( Lagos et al.|2Qlla|b ) to compute the atomic 
and molecular gas content of a series of model galaxies, and 
then used the UCL_PDR code to determine the CO luminos¬ 
ity of the molecular gas, from which they could then infer 
Aco- In their PDR calculations, it was necessary for them 
to assume some representative number density and visual 
extinction for the molecular gas. In most of their models. 


Lagos et al. (2012) assumed a characteristic number density 
n — 10"^ cm and a characteristic extinction Ay = 8. In 
these conditions, one would expect CO photodissociation to 
be completely negligible even when the ISRF is very strong, 
meaning that changes in the radiation field strength will 
only affect Aco indirectly, through its effect on the thermal 
balance of the dust. Indeed, Lagos et al. (2012) find that in 
their model, even very large changes in the radiation field 
strength have only a small influence on Aco- However, in 
practice, in realistic CMC models only a small fraction of the 
total gas mass is found in regions with a mean visual extinc¬ 
tion as high as Ay = 8 (see e.g. Clark A Glover||2014 ), and 
so by focussing only on these conditions, Lagos et al. (2012) 
ignore the large changes in the CO content and CO luminos¬ 
ity of the lower density, less shielded gas that drive much of 
the variation that we find in Aco- Lagos et al. (2012) also 
examine the effect of varying the cosmic ray ionization rate 
while keeping the radiation field strength fixed, and find that 
in this case, Aco increases slowly with increasing ionization 
rate, in general agreement with our results. 

The dependence of Aco on the UV field strength was 
also examined by Offner et al. (2014). They computed Aco 
for clouds illuminated by radiation fields with 1 and 10 

(1978) field, using the 


times the strength of the 

Draine 

3d-pdr code (Bisbas et al. 

2012). 


2012 ). When computing Aco, 


they only considered pixels in their synthetic emission maps 
with CO integrated intensities Wco > 0.45Kkms“^. They 
find that Aco = l-5Aco,gai for their Go = 1-7 run and 
Aco = 2.15Aco,gai for their Go = 17 run, a roughly 50% in¬ 
crease. This is somewhat smaller than the difference we find 
between our SFR = SFRo and SFR = 10 SFRo runs, which 
for a similar integrated intensity cut show values of Aco 
that differ by closer to a factor of two. However, it is im¬ 
portant to note that Offner et al. (2014) do not vary the 
cosmic ray ionisation rate at the same time that they vary 
the radiation field strength, and so it is likely that there is 
less destruction of CO in their Go = 17 run than in our 
SFR = 10 SFRo run. 

Finally, in a recent paper, [Bisbas, Papadopoulos A Vitij 
(2015) have studied how the typical CO abundance in a 
CMC varies as a function of the cosmic ray ionisation rate. 
They find that increasing the cosmic ray ionisation rate by 
a factor of ten or more leads to substantial destruction of 
CO and argue that this will lead to clouds in regions with 
SFR > 10 SFRo being CO-faint. However, they do not com¬ 
pute Aco for any of their clouds models and hence do not 
quantify its dependence on the SFR. 


5.3 Caveats 

There are two important methodological caveats that the 
reader should bear in mind when considering our results. 
First, as we have already discussed in Section |4.2[ we do 
not currently account for the effects of line-of-sight veloc¬ 
ity gradients when determining the effectiveness of H 2 self¬ 
shielding. This means that we will tend to overestimate the 
effectiveness of self-shielding in gas with H 2 column densi¬ 
ties of around 10 ^^ < Ah 2 < 10^^ cm“^, the regime where 
Doppler broadening of the Lyman-Werner lines dominates 
the UV absorption spectrum. Consequently, we will overes¬ 
timate the total H 2 content of the cloud. However, as most 
of the H 2 in the cloud is found in regions with much higher 
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H 2 column densities, we would not expect the omission of 
this effect to have a large impact on our results. 

Second, we have assumed that the cosmic ray ionisation 
rate is uniform throughout the cloud, or in other words that 
the cosmic ray absorption spectrum is not significantly af¬ 
fected by absorption within the cloud. Whether or not this 
is a good approximation depends on the details of the low 
energy portion of the cosmic ray energy spectrum, which 
is poorly constrained (Padovani, Galli & Glassgold 2009). 
There is some tentative evidence from astrochemical studies 
that the cosmic ray ionisation rate in dense cores may be sig¬ 


nificantly lower than in the diffuse ISM (see e.g. Indriolo k, 
McGall|2012 ), but the issue is far from being settled. If there 


is indeed a significant fall-off in the cosmic ray ionisation 
rate as we move from diffuse to dense gas, then our models 
will tend to overestimate the heating rate in the dense gas, 
and also the rate at which GO is destroyed there. The net 
effect of this on Xco is difficult to predict without detailed 
modelling. However, we note that the destruction of the rel¬ 
atively diffuse inter-clump GO component that occurs as we 
increase the SFR is a consequence of the increasing strength 
of the ISRF, and not the increase in the cosmic ray ioni¬ 
sation rate. Therefore, even if cosmic ray absorption within 
the cloud becomes significant, this behaviour will remain the 
same, and so we would expect to still see an increase in Xco 
as we increase the SFR. 

In addition to these methodological caveats, we also 
take this opportunity to remind the reader that when we 
refer to changes in the “SFR”, what we mean are changes in 
the strength of the ISRF and the size of the cosmic ray ion¬ 
isation rate, since we assume that both scale linearly with 
the star formation rate. It could well be that this assump¬ 
tion is too simplistic, and that the functional dependence 
of the radiation field strength and the cosmic ray ionisation 
rate on the local SFR is more complicated than a simple 
linear scaling. If so, then the results of our study will still 
hold, but the actual star formation rates corresponding to 
our 10 SFRo and 100 SFRo models will differ from what we 
have assumed here. 


6 CONCLUSIONS 

In this paper, we have performed a series of numerical sim¬ 
ulations of molecular cloud evolution that explore how the 
distribution of CO emission coming from the clouds changes 
as we change the strength of the interstellar radiation field 
and the size of the cosmic ray ionisation rate in the clouds. If 
one makes the reasonable assumption that these quantities 
scale linearly with the local star formation rate, then our 
results indicate how the CO-to-H 2 conversion factor for the 
clouds, Xco, depends on the star formation rate. 

We find that as we increase the radiation field strength 
and the cosmic ray ionisation rate, the CO content of our 
simulated clouds decreases. CO survives more effectively in 
well-shielded clump and filaments than in the more diffuse 
inter-clump gas, leading to a significant decrease in the filling 
factor of bright CO emission. The integrated intensity of 
the CO emission from the brightest regions increases with 
increasing SFR, owing to the heating of the gas by the higher 
cosmic ray flux, but this increase is too small to compensate 


for the loss of CO emission from elsewhere in the cloud, and 
so overall the CO luminosity of the cloud decreases. 

The change in the chemical composition of the cloud as 
we increase the SFR leads to a change in Xco - However, the 
size of this change depends on the method we use to com¬ 
pute Xco- If we consider the whole area of the cloud when 
computing the mean H 2 column density and CO integrated 
intensity, then we find that Xco changes substantially, scal¬ 
ing as roughly Xco SFR^^^ for our high-mass clouds. On 
the other hand, if we consider only regions with CO inte¬ 
grated intensities exceeding an intensity threshold IFco,cut, 
then we find that the dependence of Xco on SFR weakens 
substantially as we increase Wco,cut- The reason for this dif¬ 
ference in behaviour is that the fraction of the cloud filled 
with “dark” molecular gas (i.e. H 2 without associated CO) 
increases as we increase the SFR. This dark gas contributes 
to our calculation of Xco in the case where we do not apply 
an intensity threshold, but is ignored when we do apply a 
threshold. 

We have also explored whether increasing the density 
and turbulent velocity dispersion of our model clouds at the 
same time as we increase the SFR affects the relationship 
between Xco and the SFR. We find that if we increase the 
initial number density to no = 10 ^ cm“^ and raise the virial 
parameter to avir = 2 , then we can recover a value of Xco 
close to the canonical Galactic value even when the local 
SFR is 100 times our default value. We conclude from this 
that we will only recover a positive correlation between Xco 
and the SFR if there is not a strong correlation between 
mean cloud densities and the SFR. The fact that observa¬ 
tional determinations of Xco in starbursting systems such 
as ULIRGs find a value lower than the canonical Galactic 
value therefore provides a strong indication that the proper¬ 
ties of typical molecular clouds in these systems must differ 
significantly from those of molecular clouds in the local ISM. 
Consequently, we argue that the most important message to 
take away from this study is that it is necessary to under¬ 
stand the nature of the cloud population before one can 
properly interpret the CO emission coming from the clouds. 

Finally, we stress that if the CO emission from nearby 
galaxies is dominated by virialised clouds with masses 
around 10 ^ M© and densities around a few 100 cm“^, then 
the X-factor may vary significantly as the local SFR changes. 
In these circumstances, our standard interpretation of the 
Kennicutt-Schmidt relation, in which Xco is assumed to be 
constant throughout the galaxy, should be rethought. 
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